The Fermi-Pasta-Ulam recurrence and related phenomena for ID shallow-water waves 

in a finite basin 
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In this work, different regimes of the Fermi-Pasta-Ulam (FPU) recurrence are simulated numeri- 
cally for fully nonlinear "one-dimensional" potential water waves in a finite-depth flume between two 
vertical walls. In such systems, the FPU recurrence is closely related to the dynamics of coherent 
structures approximately corresponding to solitons of the integrable Boussinesq system. A simplest 
periodic solution of the Boussinesq model, describing a single soliton between the walls, is presented 
in an analytical form in terms of the elliptic Jacobi functions. In the numerical experiments, it 
is observed that depending on a number of solitons in the flume and their parameters, the FPU 
recurrence can occur in a simple or complicated manner, or be practically absent. For comparison, 
the nonlinear dynamics of potential water waves over nonuniform beds is simulated, with initial 
states taken in the form of several pairs of colliding solitons. With a mild-slope bed profile, a typical 
phenomenon in the course of evolution is appearance of relatively high (rogue) waves, while for 
random, relatively short-correlated bed profiles it is either appearance of tall waves, or formation of 
sharp crests at moderate-height waves. 

PACS numbers: 47.35.-i, 47.15.K-, 47.10.-g 



I. INTRODUCTION 

Nearly integrable wave systems are known to exhibit 
the Fermi-Pasta-Ulam (FPU) recurrence, when a (finite- 
size) system approximately repeats its initial state after 
some period of evolution. Starting from the first observa- 
tion of this phenomenon in the famous numerical experi- 
ment [l[ with one-dimensional (ID) lattices of nonlinear 
oscillators, the FPU recurrence and related phenomena 
were studied in many physical contexts (see, e.g., 0- 
[T3 |. and references therein). In particular, Zabusky and 
Kruskal Q discovered the solitons with their highly non- 
trivial behaviour, when numerically investigated a mech- 
anism of the recurrence for spatially periodic solutions of 
the Korteweg-de-Vries (KdV) equation. Now the theory 
of solitons has developed into one of the main branches 
of nonlinear science. 

It is a well known fact that many integrable mathe- 
matical models have their origin in the theory of water 
waves. So, the two most famous integrable equations 
are the KdV equation, first derived for weakly disper- 
sive unidirectional shallow-water waves, and the nonlin- 
ear Schroedinger equation (NLSE) which describes an en- 
velope of a train of deep-water waves fl5| . For deep-water 
waves, many analytical and experimental results concern- 
ing the FPU recurrence are known [M [H-[l|. As to 
the shallow-water regime, only some numerical studies 
for KdV and its higher-order generalizations were per- 
formed until recently (see, e.g., [IcJ El), while the FPU 
phenomenon was never considered theoretically for long 
waves in a finite flume, and also it was never studied 
experimentally in the shallow-water regime. 
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It is clear that KdV equation is not adequate for long 
waves in a finite basin where they reflect from the walls. 
Fortunately, there is another integrable model, namely 
the Boussinesq system, which approximately describes 
bidirectional shallow-water waves and therefore it is po- 
tentially useful for analytical study of the FPU recur- 
rence in a finite-length flume (concerning integrability 
of the Boussinesq system, see [20|-|23j, and concerning 
deviations of water waves from exact integrability, see 
[Tfl |24|). However, at the moment we do not have yet 
a clear theory of FPU recurrence for the shallow water, 
based on the Boussinesq system. Perhaps, a future the- 
ory should be built with the help of the sophisticated 
mathematical methods developed for obtaining spatially 
periodic solutions of integrable systems (in particular, 
for the Boussinesq model see [12|). In the present work, 
such a general purpose is not achieved, though a family 
of periodic solutions is derived here in an explicit ana- 
lytical form through a simple ansatz, which corresponds 
to a single soliton periodically moving between the walls. 
However, that solution is by no means the main result of 
our work; it just plays an auxiliary role, namely to pro- 
vide nearly "many-solitonic" initial conditions for highly 
accurate numerical experiments. 

Very recently, a short paper by the present author was 
published, where for the first time the FPU recurrence 
was studied numerically for fully nonlinear shallow-water 
waves in a finite flume [251 ] . Exact equations of motion 
for potential planar flows of a perfect fluid with a free 
surface in terms of so called conformal variables were em- 
ployed in that study [2|| [2?J . The simplest initial states 
were taken, with zero velocity field and a cosine-shaped 
free boundary. Two typical examples of recurrence are 
presented in FigsH][2J 

One of the purposes of the present work is to pro- 
vide additional numerical examples of the recurrence for 
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FIG. 1: The FPU recurrence is perfect with the initial shape 
of free surface in the form r)o(x) = 0.12ft cos(27rx/60/i): a) 
wave profiles at several time moments when the kinetic energy 
is at minimum; b) the ratio of the kinetic energy to the total 
energy; c) the maximum and minimum elevations of the free 
boundary. 



different initial states, and to demonstrate a relation of 
the FPU phenomenon in shallow-water finite basins to 
solitons of the approximate Boussinesq system. Another 
purpose is to observe what new effects appear in the dy- 
namics of long dispersive waves if the bottom boundary is 
nonuniform (it should be noted that conformal variables 
provide exact equations of motion for arbitrary nonuni- 
form bottom profile when it is parametrized by an an- 
alytical function [lEHS])- I n particular, three kinds of 
bed profiles will be considered, namely mild-slope beds, 
beds with quasi-random, relatively short-correlated cor- 
rugations, and beds with randomly placed barriers. The 
nonuniformity destroys the approximate integrability, so 
initial states in the form of several pairs of colliding 
solitons evolve to appearance of highly nonlinear wave 
events. Such steep and tall waves can be considered as 
a ID model for freak (rogue) waves sometimes arising 
in the coastal zone (the subject of freak waves is studied 
now very extensively, see [28l430l | , and references therein) . 
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FIG. 2: The FPU recurrence is less perfect with larger ini- 
tial amplitude, 770(2;) = O.Uh cos(27rx/60/i), because effects 
of non-integrability are stronger. 



II. DIFFERENT EXAMPLES OF THE FPU 
RECURRENCE 

A. Notes about numerical method 

The numerical method employed here is based on the 
parametrization of (x-periodic) free surface in terms of a 
real function p(&,t), as follows [Ml^: 

X + 1Y = Z{$ + ia(t) + {1 + iR a }p(tf, *)), (1) 

where R Q is a linear integral operator diagonal in the dis- 
crete Fourier representation, R a (m) = itanh(am). Op- 
erator R Q depends parametrically on a positive quantity 
a(t) (for details, see [2^, [27]). A fixed analytical function 
Z(Q determines conformal mapping of sufficiently wide 
horizontal stripe in the upper half-plane of an auxiliary 
complex variable C, adjacent to the real axis Im ( = 0, 
onto a region in the physical (ir, y)-plane, with the real 
axis Im £ = parametrizing the bed profile. 

Since we consider purely potential planar flows of an 
ideal fluid, the velocity field is completely determined by 
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a real function ip($,t), which is the boundary value of 
the velocity potential at the free surface. 

Exact compact expressions for the time derivatives 
pti'&it), ipti'djt), and a(t) were obtained, corresponding 
to the dynamics of potential water waves in the uniform 
gravity field a |25l27j: 



4> t = -Re[$*(T Q +i)Q] - 



(2) 

T - fl ImZ(0(3) 
(4) 



where £ = & + ia + (1 + iR a )p, $ = (1 + iR^)^ an d Q = 
(Ratp^) /\Z' (£_)^\ 2 . Operator T Q is diagonal in the dis- 
crete Fourier representation: T a (m) — — icoth(am) for 
m 7^ 0, and T a (0) = 0. Eq.© is the so called kinematic 
boundary condition at the free surface, written in terms 
of conformal variables, Eq.Q is the dynamic boundary 
condition (the Bernoulli equation), while Eq.Q takes 
into account temporal dependence of the conformal depth 
a, which is necessary for conservation of the total fluid 
volume. 

If function Z(Q) is expressed in terms of elementary 
analytical functions [such as Exp(), Log Q, and so on; for 
particular examples see (2|| H3, l3lTl33j |] , then the right- 
hand sides of Eqs.([2])-(|1]) can be easily evaluated using 
the Fast Fourier Transform routines and mathematical li- 
brary complex functions [so, in C programming language 
we have cexpO, clogO, and so on]. The above proper- 
ties form the base of the numerical method. 

In the numerical experiments, we use dimensionless 
variables (however, for graphical presentations the wave 
profiles are rescaled to a characteristic depth h = 1 m), 
and consider either flat horizontal, or 27r-periodic nonuni- 
form bed profiles [that means Z(( t +2Tr) = 2ir+Z(()} hav- 
ing an additional symmetry about the imaginary axis, 
ImZ(-(' + = ImZ(C + K"), ReZ(-C + i£") = 
— ReZ(C + i£")- This symmetry is required for simu- 
lations of waves between the vertical walls located at 
x = and at x = 7r. Of course, the functions VK^M) 
and /9($, t) should possess definite symmetries as well: 
ip(tf + 2tt, t) = t), t) = t), p{tf + 2tt, t) = 

pi'djt), p(—"3,t) = —p(i) 7 t). The symmetries are auto- 
matically preserved in time if the initial data are sym- 
metric. 

In all our simulations, the system at t = is charac- 
terized by a free surface profile y = rjo(x), and by the 
velocity field v = 0. Such initial conditions with zero 
kinetic energy E'kin were taken because they are conve- 
nient to observe the recurrence by monitoring the time 
dependence of the quantity Eki n /E, where E is the total 
energy, which is conserved in the numerical experiments 
up to 7-8 decimal digits. The function r]o(x) is even and 
periodic, therefore it satisfies the boundary conditions 
7?o(0) = rj' (L/2) — 0. A special procedure was designed, 
in order to find numerically a function p(j), 0) correspond- 
ing to a given initial profile ffa^) [25j . 
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FIG. 3: The FPU recurrence is practically absent when 
7io(x) = 0.24/i{[l + cos(27T2:/60^)] 3 /8 - 5/16}. 



B. Example when recurrence is absent 

It should be stressed that the FPU recurrence takes 
place for special initial conditions only. It is clear from 
the theoretical point of view that recurrence corresponds 
to a nearly closed trajectory on a torus in a phase space 
of an intcgrable system. The dimensionality of the torus 
is equal to the number of effectively excited degrees of 
freedom. Typically, the frequencies of that motion are 
not rationally related. Therefore, in the generic case re- 
currence is not observed. As an example when the recur- 
rence is practically absent, in Fig(3]we present evolution 
of some relevant parameters in the numerical experiment 
with T) (x) = 0.24/i{[l + cos(27ra/60fr)] 3 /8 - 5/16}. 

From this point of view, at the moment it looks some- 
what miraculous that the initial profiles in the simplest 
form f]Q(x) — Aq cos(2irx/L) demonstrate rather per- 
fect recurrences despite the number of effectively excited 
degrees of freedom can be fairly large, iV s = 5. ..7 for 
Aq/H = 0.12 and L/h = 100. ..120, as in the numerical 
experiments reported in [25| . 



C. Recurrence in the dynamics of solitons 

It is a well-known fact that intcgrable systems with pe- 
riodic boundary conditions posses so called finite-gap so- 
lutions, which are exactly finite-dimensional sub-systems 
with the dynamics on a torus. For the Boussinesq sys- 
tem, the simplest example is given in the Appendix, more 
involved cases are considered in [22| • In general case, the 
formulas are quite complicated, and it is difficult to de- 
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FIG. 4: Example of FPU recurrence when two solitons with 
different parameters are present in the system: a) wave pro- 
files at several time moments when the kinetic energy is at 
minimum; b) the ratio of the kinetic energy to the total en- 
ergy; c) the maximum and minimum elevations of the free 
boundary. 



scribe the corresponding degrees of freedom in terms of 
simple physical quantities. However, if the x-period is 
sufficiently long, an approximate description in terms of 
several colliding solitons becomes possible. In the Boussi- 
nesq model, "free" solitons are characterized by positive 
or negative dimensionless velocities s n (constant, all dif- 
ferent), and by their positions (phases) x n (t). The cor- 
responding analytical solutions are presented in [23j|. It 
is important that when two solitons with opposite veloc- 
ities si = — S2 = s collide at a position xq (or a single 
soliton collides with the wall), at some time moment the 
velocity field is identically zero along the flow domain, 
while the shape of the free surface is given by a simple 
formula T](x) — S(x — xq, s), with 



S(x — xq, s) — 



2h(s 2 - 1) 



cosh 2 [ v /3(s 2 - l)(x - x )/2h] 



(5) 



This formula was used in our numerical experiments to 
prepare initial states in the form of several pairs of col- 



FIG. 5: Two pairs of solitons with equal parameters are 
present in the system: a) wave profiles at several time mo- 
ments when the solitons are colliding; b) the ratio of the ki- 
netic energy to the total energy. 



liding solitons, placed sufficiently apart from each other. 

In a finite domain, each soliton moves between the 
walls, and additionally it acquires definite phase shift 
a(s n ) when it reflects from a wall (s n — ¥ —s n after reflec- 
tion), and phase shifts A(s„,s m ) when it collides with 
other solitons [23| . In this picture, the recurrence takes 
place when at some time moment positions of all solitons 
self-consistently return closely to their initial values. The 
simplest nontrivial example for quasi-recurrence in the 
system of two solitons is shown in FigJU 

However, the above approximate description does not 
work if we initially put several identical humps, each cor- 
responding to a pair of colliding solitons, to different po- 
sitions. FigJS] shows that the recurrence occurs in a more 
complicated way in such a case. 

If we put two or more humps, with one of them higher 
than the others, then recurrence is possible only with 
tuned values of the larger velocity. Successful examples 
are shown in FigJB] and Fig|7] 



III. EXTREME WAVES OVER NONUNIFORM 
BEDS 

We have observed through numerical experiments that 
with the flat bottom, if all the initial humps are nearly of 
the same height, then the maximum surface elevation, as 
a function of time, does not exceed much the initial value 
for a long time (not shown). What happens if the bed 
is nonuniform? This question is answered in this Section 
through a set of numerical experiments. Three typical 
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FIG. 6: Two pairs of solitons with different parameters are 
present in the system: a) wave profiles at several time mo- 
ments when the solitons are colliding; b) the ratio of the ki- 
netic energy to the total energy. 
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FIG. 7: More complicated regime of the FPU recurrence. 
Three pairs of solitons are present, with one pair stronger 
than the two other. 



examples are presented below. 

In the first nonuniform case, we simulated waves over 
a mild-slope bottom, corresponding to function Z(() in 
the form 
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FIG. 8: Rogue waves appear over a mild-slope bed, some of 
them near the wall. 



(see Fig JS]). An initial state was eight pairs of collid- 
ing solitons placed at quasi-random positions [that is, 
eight humps of the form Eq.©], with s = 1.11 (not 
shown). Over nonuniform bed, the initial solitons evolve 
to a random wave field, consisting of quasi-solitonic co- 
herent structures of different heights, together with non- 
coherent waves. After sufficiently long time, the smooth 
nonuniformity resulted in appearance of tall extreme 
wave events (see FigslSi fTUl) . when the most strong, op- 
positely propagating quasi-solitonic coherent structures 
collided. The highest wave events were observed near 
the left wall, where the depth is minimal. 

In the second nonuniform case, we took function Z(C) 
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FIG. 9: Maximum and minimum elevations of the free surface 
over the mild-slope bed. 



FIG. 11: Formation of sharp wave over a quasi-random bed. 
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FIG. 10: Big wave over the mild-slope bed near the wall: a) 
the shape of free surface, b) velocity distribution along the 
free surface. 



as a sum of several Fourier harmonics, 



Z 2 (0 = C - |jj + i £ MrC™ eMikmO, (7) 

m— 1 



with some positive coefficients C m and positive integer 
wave numbers k m . Practically, in this case the bed pro- 
file looks as a quasi-random, relatively short-correlated 
curve (see Fig llip . Again extreme waves appeared on 
a later stage of evolution, after the initial state consist- 
ing of eight pairs of colliding solitons with parameters 
s = 1.12. Comparatively to the mild slope, now ex- 
treme waves were no so high, but very sharp and typ- 
ically more asymmetric. It should be also noted that 
in this case a typical time of the transition to a ran- 
dom wave field, characterized by a rough profile of the 
free surface, was essentially shorter than for the mild- 
slope bottom. Also, quasi-solitonic coherent structures 
in the second case were relatively short-lived, while a 
non-coherent part of the wave filed was more developed. 

In the third nonuniform case, the bed inhomogeneity 
was taken in the form of 25 randomly placed nearly iden- 
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where £ n for n — 0...24 were quasi-random real numbers 
in the range from to 7r, while Cn+25 = 2n — £„ (required 
for the symmetry). The other parameters were B = 1400, 
e = 0.001, i/ = 5.0x 10~ 5 . Contrary to the previous two 
examples, the conformal mapping, corresponding to func- 
tion ^(C), has singularities in the upper half-plane of the 
complex variable C; however, all the singularities are far 
enough above the free surface. The large value of the 
parameter B allowed us to take into account only several 
nearest barriers when the function Z{Q was numerically 
evaluated, thus significantly reducing the computational 
cost comparatively to many different possible choices for 
barrier shapes. 

Eight initial pairs of solitons had s = 1.11. In this 
numerical experiment, extreme waves appeared as well 
(see Fig fT^j) . Some of them were quite tall, while other 
had a moderate height, but a very sharp crest (these 
were essentially asymmetric). In general, the third case 
is much similar to the second case (compare Fig |12l and 
FigHU. 
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FIG. 12: Formation of extreme waves over a bed with ran- 
domly placed barriers. 



IV. SUMMARY AND DISCUSSION 

In this work, it has been demonstrated through highly- 
accurate numerical simulations of exact equations of mo- 
tion for planar potential flows of a perfect fluid with 
a free surface, that shallow-water dispersive waves with 
moderate amplitudes A/h < 0.12 can exhibit the Fermi- 
Pasta-Ulam recurrence in a finite basin for various ini- 
tial states. However, the best quality of the recur- 
rence is observed for the initial free surface in the form 
t]q(x) — Aq cos(2ttx/L) 1 and ipo(x) — 0. In that spe- 
cial case, velocities of all arising solitons appear self- 
consistently tuned to their phase shifts in mutual col- 
lisions, which property results in remarkably perfect re- 
currence to the initial state even in quite long flumes. A 
mathematical reason for such self-consistency is not clear 
at the moment. The FPU quasi- recurrence is also robust 
with initial states in the form of two solitons. For larger 
number of solitons, quasi-recurrence is possible with spe- 
cial values of parameters only. 

All our numerical results are based on the inviscid the- 
ory. Of course, in the reality a viscous friction will act 
against the recurrence. However, in [25| it was estimated 
that a relative effect of the viscous friction near the bot- 
tom and near the side walls of the flume becomes small 
if all the spatial scales increase proportionally. 

In the quasi-integrable regime over a flat horizontal 
bed, with initial states in the form of several nearly equal 
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solitons, formation of extreme waves appears effectively 
suppressed, since the solitons preserve their strengths for 
a long time. When the approximate integrability de- 
stroyed by a bed nommiformity, the system evolves to 
a random-wave-field regime where quasi-solitonic coher- 
ent structures of different amplitudes are present, some of 
them being stronger than the initial solitons. When the 
strongest oppositely propagating structures collide, fairly 
extreme waves arise. The most high extreme waves were 
observed for a mild-slope bed profile, while for relatively 
short-correlated bed inhomogeneities the extreme waves 
were typically less tall but more sharp-crested. Similar 
effects were observed both for waves between the vertical 
walls, and for waves with periodic boundary conditions 
without the additional symmetry (not shown). 

Our present results for extreme events in bidirectional 
wave fields over nonuniform beds may have some rele- 
vance to the problem of rogue (freak) waves in coastal 
zone, but only if the coast is in the form of a wave- 
reflecting cliff rather than a wave-absorbing beach. 
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tion for Support of Young Scientists and Leading Scien- 
tific Schools (project no. NSh-6885.20f 0.2), and by the 
Presidium of the Russian Academy of Sciences (program 
"Fundamental Problems of Nonlinear Dynamics"). 



Appendix: A single soliton between the walls in the 
Boussinesq model 

The Boussinesq equations for weakly nonlinear, weakly 
dispersive long water waves in the dimensionless variables 
[rj/h —> 7], \[3xj2h — > x, \j3gjhtj2 — > t] take the form 



u t + uu x + r\ x = 0, 

■qt + [(1 + T))u] x + -u xxx 



0. 



(A.f) 
(A.2) 



where rj is the vertical displacement of the free surface, 
and u ==_0x is a horizontal velocity. Following Zhang 
and Li [23j, we transform the above system to a more 
symmetric form 



qt + T^qxx + q 2 r = 0, 



-U + -r xx + r q = 0, 



(A.3) 
(A.4) 



where new real unknown functions q(x,t) and r(x,t) ex- 
press the old functions rj(x,t) and u(x,t) in the following 
manner: 



u = q x /q, V = -f + qr + u x /2. 



(A.5) 



We note the system of equations (|A.3|) and (|A.4[) is for- 
mally similar to the focusing nonlinear Shroedinger equa- 
tion 2iip t + i> xx + 2ip 2 i/j* = and its complex conju- 
gate — 2iij)t + i/)* x + 2tp* 2 i[> = 0. Therefore we can ap- 



ply a simple generalization of the Akhmediev-Eleonskii- 
Korneev-Kulagin ansatz || [U, [35[ , and search for a so- 
lution of Eqs. (fA~3)) -(|A~4 | in the form (q,r) = [U(x,t) ± 
y/Z(t)] exp(±P(i)). Indeed, by doing so one can obtain 
and integrate a system of equations for the unknown 
functions U(x,t), Z(t), and P(t). At some point, the 
problem is reduced to the analysis of two equations (com- 
pare to dj]): 

Z 2 - 16Z 4 + 16wZ 3 - 4(ft + w 2 )Z 2 - AbZ = 0, (A.6) 



U 2 + U 4 + 2(w ~3Z)U 2 + 2-^=U + {2wZ - 3Z 2 - b) = 0, 

v Z 



(A.7) 

where w, h, b are some constants (there is also the third 
equation P + 2Z = w). 

However, we prefer not to deal with a function of two 
variables as U(x,t), and therefore we use here a slightly 
less general ansatz which still contains physically inter- 
esting solutions, and where the variables are separated 
from the very beginning: 

Q(t) 



q{x,t) = F(t) + 
r(x,t) = G(t) + 



D{x)+A{ty 
D{x) + A{t) ' 



(A.8) 
(A.9) 



We take the only x-dependent function D(x) satisfying 
the relations 



D 2 X = ^ 2 (D 2 -l)(l-e 2 D 2 ) 
= cD 2 -5D A -P, 
D TT = cD — 2SD 3 . 



(A.10) 
(ATI) 

Thus, it is one of the elliptic Jacobi functions (for their 
definitions and properties, see, e.g., [36j]): 



D(x) = nd(2(ix, sj\ - e 2 ). (A.12) 
The x-period of this function is L = 1(e)/ p,, where 

dz 



(A.I3) 



n ^(z 2 -l)(l-€ 2 z 2 ) 

Now we substitute the ansatz (|A.8I) - (|A.9[) into the 
system (|A.3|) - (IA.4I) . Using relations (IA.10[) and (|A.I1I) . 

we obtain the following set of equations (which are co- 
efficients in front of different powers (D + A)~ n , for 
n = 0, 1, 2, 3, or their linear combinations): 







-6AQ = 


0, 


(A.14) 




-G + FG 2 ■+ 


- SAR = 


0, 


(A.I5) 


-66 A 2 ) + 2QFG4 


RF 2 = 


0, 


(A.I6) 


*+!(. 


- 66 A 2 ) + 2RFG 4 


-QG 2 = 


o, 


(A.I7) 




2A + QG 


- RF = 


0, 


(A.18) 




(cA~26A 3 ) + QG 


+ RF = 


0, 


(A.I9) 




(cA 2 - 6 A 4 - /3) 


+ Qi? = 


0. 


(A.20) 
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(The last equation actually appears twice.) It is easy 
to show that the two algebraic relations are consistent 
with the five differential equations. It follows from these 
equations also that 



FG = 5 A 2 + 7, 



(A.21) 



where 7 is a constant. Now we take the squared equation 
for A and obtain 4i 2 = (RF - QG) 2 = (RF + QG) 2 - 
AFGQR = (cA - 2SA 3 ) 2 + 4(7 + SA 2 )(cA 2 - SA 4 - /3), 
that is an easily solvable first-order equation 



4i 2 



-47M 4 + A 2 (c 2 + 4c7 - 4/36) - 4 7 /3 



4 7 <5(A 2 -a 2 )(a 2 -A 2 ), 



(A.22) 



where 



"1,2 



1 



C 
47(5 




4jS 



8' 
(A.23) 

The solution of Eq. (|A~22|) is again expressed through an 
elliptic function: 



Alt) = a% nd ta-i y jS, \ 1 — 



«2 



n 2 



= «i nd (£, k) . 



(A.24) 

Since, by definition, nd(£, k) = l/dn(£, k) (see |3g|), and 
[dn(£, = — K 2 sn(£, k) cn(£, «), for the time derivative 
A(t) we have 



(A.25) 



A(t) = u\U2\J "f5n sd(£, k) cd(£, k). 



Thus, we have obtained explicit expressions for the 
quantities A, QR, FG, and (RF + QG). From these 
one can also extract the relations Q/F and R/G, since 



R Q _ {RF + QG) _ {cA - 2SA 3 ) 
G + F~ FG ~ SA 2 + 7 
R Q _ 2A 2i 
G ~ F ~ FG ~ <L4 2 + 7 ' 



, (A.26) 



(A.27) 



The obtained information is sufficient to construct the 
velocity u(x, t) and the free surface elevation rj(x,t) via 
formulas (|A.5j) . because 



1 



1 



(D + A + Q/F) (D + A) 
-l + FG> {RF + QG) > 
-u x /2. 



D x , (A.28) 
QR 



(D + A) (D + Af 



(A.29) 



Using relations c = 4/j 2 (1 + e 2 ), 5 = 4^i 2 e , and f3 = 4// , 
it is possible to show that a.\ < 1 and cti > 1/e, and 
therefore at definite time moments function A(t) takes 
values A\ — 1 or A2 = Simultaneously, at those 

time moments either Q = 0, or R = 0. When Q = 0, 
then the velocity field u(x) is zero everywhere, while the 
free surface profile is either r)x(x), or T)2(x), where 



■qi (x) = -1 + 7 + 4/j 2 
m(x) = r] 1 (x-L/2), 



(1-e 2 ) 



D(x) + 1 



, (A.30) 
(A.31) 



We see that the best choice for the constant 7 is 7 = 
1 — 4/j 2 e, since in this case 7?i m i n = 0. Function 771(2:) has 
a single hump at x = 0, and thus it corresponds to the 
moments when a soliton is colliding with the left wall, 
while 772(2;) corresponds to the collisions of the soliton 
with the right wall at x = L/2. 

In the limit t « 1 we have L — > 00, and D(x) w 
cosh(2/jx), so 



4^ 2 



2^ 2 



1 + cosh(2/ix) cosh 2 (/ia;) 



(A.32) 



The full solution in this limit is given by the following 
formulas (it is interesting to note the solution below is 
essentially Eq.(52) of Ref.[34j for the focusing nonlinear 
Schroedinger equation, but evaluated at imaginary time): 



q(x,t) = 
r(x,t) = 



2^ 2 cosh(2^ v / 1 + M 2 ) - 2^1 + y 2 sinh(2^ vT + fi 2 ) 
y/l + y 2 cosh(2^a;) + cosh(2/j/jyT+/?) 



2fjL 2 cosh(2t/J V / l + A* 2 ) + Z^V 1 + [i 2 sinh(2f/x >/l + /i 2 ) 
yl + /i 2 cosh(2/iir) + cosh(2rj/xyl + /i 2 ) 



(A.33) 
(A.34) 



Collision of the soliton with the wall at x = occurs at It is easy to derive that before and after the collision, 
t = t*, when the soliton (at x > 0) moves with the velocities s = 

T\A + M 2 , respectively. 



sh(2<^ v /l + M 2 ) = \j\ + /i 2 . 
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